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1.  Introduction 

As  zeros  of  analytic  functions,  branch  points  are  a  well  known  and  much  studied  in  complex 
analysis  [1].  The  character  zero  in  amplitude  is  accompanied  by  a  2n  circulation  in  the  phase 
and  a  choice  of  the  principal  domain  which  establishes  the  branch  cut.  It  has  been  found  that 
under  the  right  conditions  electromagnetic  fields  contain  branch  points,  a  monochromatic  wave 
originating  as  a  point  source  then  propagating  through  atmospheric  turbulence  being  one  such 
example.  It  has  been  found  that  branch  points  pose  a  difficult  problem  for  adaptive  optics  [2], 

While  the  necessary  and  sufficient  conditions  for  branch  point  formation  are  not  known, 
generally  branch  points  begin  to  form  when  the  (5 /6)th  distance  (propagation)  moment  of  the 
atmospheric  structure  function  exceeds  0.1  [3].  Hence,  the  distance  to  the  onset  of  branch  point 
formation  varies  depending  on  the  atmospheric  strength.  The  propagation  is  important;  without 
it,  atmospheric  branch  points  would  not  form.  Here,  we  are  interested  in  how  the  propagation 
affects  the  creation  and  evolution  of  branch  points. 

Because  we  are  interested  in  astronomical  imaging,  we  restrict  ourselves  to  the  paraxial 
propagation  region  with  the  axis  defined  by  the  optic  axis  of  the  astronomical  telescope.  The 
main  mirror  of  the  telescope  defines  a  plane  orthogonal  to  the  direction  of  propagation;  we  call 
this  ‘the  pupil  plane’.  A  wave  is  said  to  be  ‘in  the  pupil’  if  it  is  evaluated  at  this  plane.  The 
pupil  plane  is  the  area  of  interest  for  the  astronomical  community  since  this  is  the  part  of  the 
wave  that  forms  the  image.  While  ultimately  we  are  solely  interested  in  the  wave  in  the  pupil, 
the  intervening  region  between  the  turbulence  and  the  pupil  is  the  main  concern  here.  Also,  as 
is  well  known,  the  propagtion  the  propagation  is  described  by  the  Fresnel  integral. 

By  the  nature  of  the  telescope  in  astronomical  imaging,  the  optic  axis  is  aligned  with  the 
direction  of  propagation.  As  such,  we  define  %  to  be  the  direction  of  propagation.  The  pupil 
becomes  a  two  dimensional  x  —  y  slice  of  the  three  dimensional  propagating  wave.  Later  in 
this  study,  we  will  consider  a  sequence  of  these  transverse  planes  along  the  z  axis  as  the  wave 
propagates  from  the  turbulent  region  to  the  pupil. 

Once  in  the  pupil,  the  wave  is  conditioned  using  adaptive  optics  (AO).  A  conventional  AO 
sytem  only  correct  the  phase.  It  does  so  by  subdividing  the  telescope’s  pupil  into  a  hexagonal 
array  of  subapertures.  Within  each  subaperture,  a  phase  gradient  measurement  is  made  via  a 
Shack-Hartmann  wavefront  sensor,  and  a  deformable  mirror  actuator  is  placed  at  the  corner  of 
each  subaperture.  The  gradient  measurements  are  stitched  together  using  a  least  square  recon¬ 
structor.  It  has  been  shown  that  the  least  square  reconstructor  does  not  measure  branch  points; 
however,  other  reconstructors  do  [4]  and  we  assume  that  one  of  these  is  used. 

Previous  studies  in  conventional  adaptive  optics  solely  considered  branch  points  in  the  pupil 
plane,  a  single  transverse  slice  of  the  propagating  wave.  Such  a  restriction  comes  about  naturally 
from  the  nature  of  the  wavefront  sensor  measurement;  by  its  construction  a  wavefront  sensor 
measures  in  a  single  transverse  slice,  always  conjugate  to  the  deformable  mirror  and  sometimes 
conjugate  to  the  pupil.  In  this  frame,  the  branch  points  are  assumed  to  be  independent,  although 
after  the  frame  is  reconstructed,  branch  cuts  are  placed  between  branch  points  typically  using 
a  minimum  distance  metric  [5].  Here,  we  consider  propagation  prior  the  reaching  the  pupil  and 
enumerate  on  how  this  formation  impacts  the  pupil  plane  wavefront  sensor  measurement. 
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The  ultimate  goal  of  this  study  is  to  extend  the  state-of-the-art  of  adaptive  optics  such  that  it 
can  compensate  for  atmospheric  turbulence  even  when  branch  points  are  present.  In  this  paper 
we  demonstrate  a  narrow  but  crucial  result:  atmospherically  created  branch  points  must  appear 
in  pairs  of  opposite  polarity  and  at  the  instant  they  appear,  they  have  infinitesimal  separation. 
In  the  later  papers  in  the  series,  we  will  expand  on  this  results  with  the  eventual  demonstration 
that  a  distribution  of  branch  points  measured  in  the  pupil  plane  of  a  telescope  can  be  used  to 
estimate  the  three  dimensional  turbulence  that  created  that  distribution. 

Although  this  initial  demonstration  is  done  under  the  assumptions  of  Kolmogorov  turbulence 
and  Taylor’s  frozen  flow  hypothesis,  these  restrictions  are  not  severe.  Assuming  Kolmogorov 
turbulence  is  a  weak  assumption  since  almost  all  atmospheric  turbulence  is  described  by  it  [6], 
Assuming  frozen  flow  merely  places  an  upper  bound  on  sampling  rate  since  any  continuous 
function  will  appear  linear  when  sampled  fast  enough.  Later  papers  will  relax  this  assumption 
and  allow  for  evolving  turbulence.  In  passing,  it  is  important  to  note  that  we  are  not  interested 
in  lasers  or  cavity  generated  waves.  Also,  since  we  assume  that  the  wave  originates  as  a  point 
source,  speckle  is  also  not  an  issue. 

To  this  end,  after  a  short  background,  we  present  in  Section  3  the  basis  of  the  problem. 
In  Section  4,  we  demonstrate  a  means  to  mathematically  capture  the  notion  of  causality  in  the 
propagation  apropos  to  creation  of  branch  points.  Following  which.  Section  5  demonstrates  that 
upstream  of  the  pupil,  atmospherically  created  branch  points  are  created  in  pairs  infinitesimally 
close  together  and  having  opposite  polarity.  The  latter  part  of  Section  5  discusses  both  propa¬ 
gation  of  creation  pairs  and  creation  of  large  numbers  of  branch  points  in  the  pupil.  Sections  6 
and  7  discuss  how  branch  cuts  evolve  in  distance  and  how  choice  of  the  branch  cut  affects  the 
principal  domain.  Next,  in  Section  8,  we  speculate  on  the  nature  of  the  proof  to  include  that  the 
branch  points  must  have  both  the  velocity  and  distance  to  the  turbulence  layer  encoded  in  their 
statistics.  Finally  Section  9  wraps  up  the  work. 

2.  Background 

There  has  been  effort  in  the  adaptive  optics  community  both  to  measure  branch  points  and  to  un¬ 
wrap  the  phase  in  the  presence  of  branch  points  [2,4,7-12].  Unwrapping  the  phase  is  of  funda¬ 
mental  importance  because  AO  wavefront  sensors  eith  measure  gradient  (e.g.  Shack-Hartmann) 
or  modulo  2k  phase  (e.g.  Self-referenced  Interferometer)  the  deformable  mirrors  which  correct 
the  atmospheric  disturbance  are  continuous  face  sheet  devices  and  can  only  correct  continuous 
phase.  Even  with  this  unwrapping.  Fried  and  Vaughn  [2]  demonstrated  that  there  is  ambiguity 
in  assigning  a  continuous  phase  function  with  a  scalar  field  due  to  the  presence  of  branch  points. 
These  works  make  the  implicit  assumption  that  branch  points  are  independent  entities,  i.e.  they 
appear  singely  and  can  be  treated  as  such. 

In  his  seminal  paper.  Fried  [7]  demonstrated  that  the  gradient  of  the  phase  is  composed  of 
the  gradient  of  a  smooth  function  (the  scalar  potential)  and  the  curl  of  a  discontinuous  function 
(the  vector  potential)  shown  to  be  the  Hertz  potential.  There  has  been  work  following  Fried’s 
foundational  work  [4, 8, 9, 1 1]  and  one  [10]  which  demonstrated  that  the  direction  of  the  phase 
map  (gradient)  plays  a  role  in  the  setting  of  branch  points.  For  most  of  AO  history,  circulation 
in  the  phase,  formerly  called  ‘slope  discrepancy’,  was  caused  only  by  noise.  This  is  because 
in  near-zenith  imaging,  the  phase  is  irrotational,  that  is,  there  are  no  branch  points.  Brennan 
[13]  demonstrated  that  slope  space  has  as  a  basis  the  slope  discrepancy  subspace  and  the  least 
squares  subspace.  In  the  regime  in  which  we  work,  slope  discrepancy  consists  of  both  noise 
and  branch  points. 
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3.  Problem  Formulation 


Assume  there  is  a  monochromatic  point  source  beyond  the  atmosphere  such  that  z  <  0.  Prop¬ 
agate  it  as  an  expanding  spherical  wave  until  it  reaches  the  atmosphere  at  zatm,  i.e.  0  <  za tm  < 
Zpupii  where  it  encounters  the  atmosphere.  The  finite  extent  of  the  atmospheric  acts  as  an  apodiz- 
ing  element  giving  the  wave  finite  support,  Clatm.  Immediately  prior  to  satm,  the  wave  has  con¬ 
stant  amplitude  and  a  spherical  phase,  (j>sph(za tm)-  At  the  atmospheric  layer,  as  the  beam  prop¬ 
agates  through  the  atmosphere’s  varying  index  of  refraction,  n(ralm),  the  phase  of  the  beam 
takes  phase  kratm  with  A  the  wavelength  of  the  light,  and  ratm  the  coordinates  transverse  to  the 
direction  of  propagation.  The  phase  in  the  plane  zatm  is  then 

2n  , 

0Zatm  —  n\ratm )  +  tyslip 

yielding  for  th  wave  at  zatm 

f{?atm ■  Zatm )  =  (1) 


with  A( rafm,zatm)  constant.  The  wave  propagates  from  zatm  until  it  reaches  the  pupil  at  zpupii 
where  it  is  measured. 

The  wave  for  any  plane  z  >  ’aim  is  given  by  the  Fresnel  propagator. 


/(r.z) 


(2) 


where  A  the  wavelength  of  the  light,  k  =  ^  the  wavenumber,  £latm  the  aperture  through 
which  the  wave  passes  at  plane  zatm,  z  the  coordinate  location  along  the  direction  of  propa¬ 
gation,  r  =  (x,y)  are  the  coordinates  transverse  to  the  direction  of  propagation,  ratm  the  trans¬ 
verse  coordinates  at  the  atmospheric  interaction,  /( r,  0)  the  field  at  the  turbulence  layer,  and 
/(r.z)  the  field  a  distance  z  downrange  from  the  turbulence  layer.  As  is  well  known  from 
the  Plancherel  Theorem  and  the  Riemann-Lebesgue  Lemma,  under  these  conditions  at  each 
Z,  /(r,z)  €  C°°(R2)  almost  everywhere  [14],  i.e.  /  is  an  element  of  the  set  of  continuously 
differentiable  functions  on  R2  except  on  sets  of  measure  zero. 

As  the  wave  continues  to  propagate,  it  reaches  the  pupil  at  z  =  zpupii  and  is  measured  by  an 
adaptive  optics  wavefront  sensor  on  a  telescope.  This  measurement  is  as  discussed  in  [7]  and 
the  associated  difficulties  in  estimating  <p-plIllll  from  these  measurements  as  discussing  in  [2].  At 
ZpUpii,  the  telescope  both  apertures  the  wave,  applying  the  finite  support  ff-pupM ,  and  also  removes 
the  spherical  component,  ^sp/Zzpupii),  of  the  wave.  The  wave  can  then  be  written 


/(t*? Zpupii)  —  ^ZpupilA(rjZpupil)e  ^ PuPl1 


(3) 


with  A (r, Zpupii)  non-constant,  0-pllpil  the  phase  in  the  pupil  after  removal  of  the  spherical  focus 
term  by  the  telescope,  and 


Q. 


Zpupii 


1  r  <D 

0  r  >D 


where  D  is  the  diameter  of  the  telescope.  Our  interest  here  is  in  determining  how  cTpupi  came  to 
be,  and  this  is  determined  by  the  intervening  propagation. 

So  following  Fried  [7, 1 1],  in  any  transverse  plane  along  z  with  zatm  <  z  <  zpUpii.  let  w(x,y)  := 
/(r,z)  be  the  scalar  function  representing  an  optical  field  in  plane  z.  A  scalar  function  can  be 
assumed  because,  although  it  makes  calculation  easier,  it  does  not  affect  the  results.  Let  u(x,y) 
and  v(x,y)  be  the  real  and  imaginary  parts  of  w(A.y),  respectively,  with 


w{x,y)  =  u(x,y)  +  iv(x,y). 


(4) 
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The  amplitude,  A  (x,y),  and  phase,  fi(x,y),  of  this  field  are  given  by 


(5) 

(6) 


A(x,y)  =  ^Ju2(x,y)  +  v2(x,y), 
<Kx,y)=arg(w(x,y))+2nK, 


where  on  physical  grounds  0  is  a  continuous  function  [11]  and  this  implies  K  =  K(x,y)  with  K 
consisting  entirely  of  two-dimensional  step  functions  of  height  2n.  Since  w(x,y)  :=  f(r,z)  € 
C°°(R2)  almost  everywhere,  w(x,y)  is  continuously  differentiable  over  regions  that  do  not  in¬ 
clude  the  discontinuities  inherently  defined  by  K.  In  this  way,  we  can  write 

w(x,y)  =  A[x,yy^  e  C2(K2\B)  (7) 

with  the  notation  C2(Mr\B)  denoting  the  set  of  continuously  differentiable  functions  on  R2 
except  on  the  set  of  points  given  by  set  B  which  is  the  set  of  coordinates  along  which  the  branch 
cuts  reside.  Both  A(x,y)  and  fi(x.y)  are  continuously  differentiable  on  C2(M.2\B). 

It  was  shown  by  Fried  [7]  that  phase  found  in  adaptive  optics  caused  by  atmospheric  turbu¬ 
lence  is  given  by  is  a  smooth  component  and  a  discontinuous  component  due  to  branch  points 
and  can  be  written 

§  pupil  —  films  T  filiid  (8) 

where  fipupi •/  is  the  phase  of  the  wave  in  the  pupil,  fii„ts  is  the  smooth  component,  and  fi^  is  the 
discontinuous  component  due  to  branch  points.  The  branch  point  component  is  given  by 


0hid(r)  =  Im 


nLi  (x~xk)  +  i(y-yk)  1 
nf=1(x-x'i)+i(y-y'l)  \  J 


(9) 


with  (xicyit)  the  coordinates  of  branch  points  with  positive  polarity,  (x,l,y'!)  the  coordinates  of 
branch  points  with  negative  polarity,  K  the  number  of  branch  points  with  positive  polarity,  and 
K'  the  number  of  branch  points  with  negative  polarity. 

From  the  theory  of  analytic  functions,  a  branch  point  at  (xp,yp)  is  equivalent  to  both  a  zero 
in  the  amplitude,  w(xp,yp )  =  0,  and  a  circulation  in  phase.  The  circulation  in  phase  is  given  by 


JO 


dl'  •  =  ±2n 


(10) 


for  C'  a  closed  curve  encircling  the  one  branch  point  at  (xp.yp).  and  V0  the  gradient  operator 
acting  on  fi.  There  are  many  ways  of  writing  Equation  10,  but  the  form  chosen  here  mimics 
how  the  branch  points  are  measured  by  an  AO  system’s  wavefront  sensor.  Also,  for  discussional 
clarity,  we  call  sign  ^  fc,  dV  •  V ^(x1  ,y')  =  ±2^^  the  polarity  of  the  branch  point.  Also,  it  is  well 
known  that  the  zeros  of  an  analytic  function  are  isolated  and  so  the  set  of  all  branch  points  in 
w(x,y)  are  also  isolated.  It  could  be  said,  in  a  one  sentence  description,  that  this  paper  is  about 
how  zeros  of  the  analytic  wave  are  created  and  evolve. 

As  is  well  known  [1],  mappings  with  branch  points  are  multi-valued  and  in  order  to  create 
a  single  valued  function,  a  principal  domain  must  be  chosen.  Choosing  a  principal  domain 
establishes  the  location  of  the  branch  cut,  and  hence  the  set  B.  Along  the  branch  cut,  a  2n 
discontinuity  exists.  Typically  in  complex  analysis,  the  branch  cut  (and  hence  B )  can  be  chosen 
arbitrarily,  and  when  calculating  residues,  it  is  chosen  to  minimize  the  difficulty  in  integration. 
Here,  we  will  show  that  this  cannot  be  done,  i.e.  that  the  location  of  the  branch  cut  (and  hence 
the  principal  domain)  is  determined  by  the  physical  system;  this  will  be  clarified  in  Section  7. 

To  establish  notation,  for  branch  point  p,  let  Bp  =  {(x,y)c}cGbranchcutofP  be  the  set  of  coor¬ 
dinate  values  along  which  the  branch  cut  lies.  Define  B(x,y;z)  to  mean  “at  each  propagation 
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distance  z,  B(x,y;z )  =  {(x,y)c}cgaii  branch  cuts  is  the  set  of  coordinates  of  all  branch  cuts  in  the 
plane  transverse  to  z”.  For  a  wave  propagating  in  z,  no  apriori  restrictions  are  placed  on  varya- 
tions  in  B  although  it  will  be  demonstrated  later  that  B  is  greatly  constrained  by  causality. 

Finally,  in  order  to  study  the  regions  where  branch  points  form,  we  use  uniform  convergence 
of  a  sequence  of  functions.  Uniform  convergence  is  covered  in  standard  textbooks,  for  instance 
Rudin  [15].  Heuristically,  uniform  convergence  occurs  when  a  sequence  of  functions,  /„,  con¬ 
verges  to  a  function  simultaneously  at  all  coordinate  values.  Precisely,  given  a  sequence  of 
functions  {/«}»=i,2,3,...  on  a  set  R2.  {/„}«= 1,2,3,...  converges  uniformly  on  R2  if  for  every  e  >  0, 
3N  s.t.  n  >  N  =>  |  fn(x)  —  f(x)  \  <  e,  Vx  G  R2.  Here  for  ease  of  discussion: 

A  wave  is  said  to  propagate  uniformly  if  and  only  if  {/„}  — >  /  uniformly. 

Uniform  propagation  is  equivalent  to  enforcing  causality,  i.e.  to  preventing  a  wave  from  chang¬ 
ing  instantaneously  across  all  space  when  propagating  from  cross  sectional  plane  z  to  cross 
sectional  plane  z  +  £  where  e  is  an  arbitrarily  small  distance.  Interestingly,  we  will  show  that 
after  interaction  with  a  turbulence  layer,  propagation  is  uniform,  and  this  by  itself  is  sufficient 
to  make  statements  about  how  branch  points  appear  in  the  pupil. 

4.  Fresnel  Propagation  is  Uniform 

Say  there  is  a  monochromatic  wave  propagating  along  £  and  say  the  telescope  pupil  exists 
at  plane  zpUpii  -  We  study  the  behavior  of  waves  as  they  move  forward  in  z.  So,  create  a  se¬ 
quence  in  z,  labeled  z„,  such  that  { z„ }  is  monotonically  increasing  and  z,n  —■ ►  zpupii  as  11  — >  °°. 
Let  w(x,y)  =  /(x,y,zpupii)  be  the  two  dimensional  scalar  function  representing  a  cross  section 
of  the  traveling  wave  in  the  pupil.  Then  create  a  set  of  two  dimensional  functions  of  transverse 
planes  approaching  plane  zpupii  in  steps  of  2  by  using  f(x,y,zn)  with 

Zn  • —  £pupil  U  G  7L,  (1 1) 

n 

and  to  preclude  the  un-physical  occurrence  of  propagating  back  in  time,  assume  without  loss 
of  generality  that  z  >  1.  Note  in  passing  that  the  choice  of  {z,n }  is  arbitrary  as  long  as  zn  is 
monotonically  increasing  with  n  and  z„  — >  zpupii-  Then,  let 


fn  ■ =f(x,y,Zn ) 


(12) 


and  create  a  set  {/„}. 

First,  assume  that  the  first  branch  points  form  after  the  pupil,  i.e.  for  z  <  zpupii  there  are  no 
branch  points.  (This  restriction  will  be  lifted  in  a  subsequent  Section  5.)  Hence,  0/„v/  =  0  and 

{0AO  =  $[ms 

(13) 

B(x,y;z<  zpupii)  =0 


with  0  denoting  the  null  set. 

It  must  be  shown  that  the  sequence  {/„}  is  such  that  /„  — >  /(r,z)  uniformly.  Using  results 
from  functional  analysis,  this  is  easy,  if  opaque.  Following  [14,  Ch.  IX],  given  Jz?(R")  the 
Schwartz  space  of  C°°  functions  of  rapid  decrease,  Jz ?(R”)  is  a  Frechet  space  and  hence  is 
complete,  and  further  the  Fourier  transform  is  a  linear  bicontinuous  bijection  from  Jz ?(R") 
onto  Since  cos()  €  C°°(R”)  and  sin()  G  C°°(R”)  with  1 1  •  |  |o»  <  1,  given  a  function, 

g  €  j£f(R"),  multiplication  by  either  the  sine  or  cosine  also  yields  a  function  in  «5f(R").  Hence, 
the  Fresnel  transform  given  by  Equation  2  is  also  a  linear  bicontinuous  bijection  from  Jz? (R2) 
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onto  Jf(M2).  So,  given  a  function,  /  G  Jf(K2),  there  exists  a  sequence  {/„}  with  /„  G  Jz? (K2) 
such  that  /„  — >  /.  Since  the  Fresnel  transform  is  linear,  choosing  /„  as  in  Equation  12  with  z„  as 
in  Equation  11  yields  the  desired  result,  /  — >  /„  uniformly  [14],  A  more  lucid  demonstration, 
which  only  uses  properties  of  the  Fresnel  integral,  is  given  in  Appendix  A. 

This  shows  that  when  no  branch  points  are  present,  f,  — »  f(r,z)  uniformly  on  R2\0.  More¬ 
over,  since  analytic  functions  are  unique  and  given  Equation  7,  this  implies  that  both  the  am¬ 
plitude  and  phase  propagate  uniformly,  i.e.  create  sets  {A,,}  and  {(j),,}  with  An  =  A(x,y,zn) 
and  (j)n  =  0(x,y,z„ )  and  z,n  =  zpupii  —  ^  as  before.  Then  given  that  fn  — >  /(r,z)  uniformly , 
A„(r)  — >  A(r)  uniformly  and  </>„( r)  — >  <j>(r)  uniformly. 

This  is  the  crux  of  the  matter;  the  rest  that  is  presented  here  follows  easily  from  this  result. 
Interestingly,  this  by  itself  is  sufficient  to  make  statements  about  how  branch  points  appear  in 
the  pupil;  how  this  is  so  is  shown  in  the  subsequent  sections. 


5.  Paired  Branch  Points  and  Separation 

To  prove  that  branch  points  are  only  created  infinitesimally  close  together  in  pairs  of  opposite 
polarity,  first,  in  Section  5.2,  it  is  shown  that  pairs  of  branch  points  of  opposite  polarity  can  be 
created  without  violating  causality  but  only  when  created  infinitesimally  close  together,  while 
in  Section  5.3  demonstrates  this  is  false  for  odd  numbers  of  branch  points. 


5.1.  Arbitrary  Non-trivial  Branch  Point  Phase 


We  will  demonstrate  that  paired  branch  points  can  occur,  but  only  when  created  infinitesimally 
close  together. 

So,  suppose  we  have  a  wave  propagating  in  the  z  direction  such  that,  prior  to  the  plane 
zpupii,  the  wave  does  not  contain  branch  points.  Suppose  further  that  in  plane  zpupii  an  arbitrary 
distribution  of  branch  points  appear,  i.e.  for  z  <  zpupii-  fitud  =  0  and  for  z  >  zpupii,  <l>i,id  7^  0  for 
which 

{•fepupil  =  Vims 

(14) 

B(x,y;z  <  zpupit )  7^  {} 

and 

{•fepupil  =  Qlms  +  Qhid 

(15) 

B(x,y;z  <  zpupii)  7^  {} 

Let  f(x,y,z)  denote  the  wave  prior  to  plane  z  and  create  the  sequence  {/„}  with  /„  = 
f(x,y,zn )  as  before  by  using  the  sequence  z„  =  zpupii  —  jj.  Let  /(x,y,zpupii)  denote  the  wave 
in  the  pupil.  Let  w(x,y)  =  f(x,y,zpup ii)  be  the  two  dimensional  scalar  function  representing  a 
cross  section  of  the  traveling  wave  in  the  pupil,  and  as  was  shown  by  Fried, 

w(x,y)  =A(x,)0e,(*“(**)+*I"(*J')  (16) 


For  the  general  case  of  an  arbitrary  distribution  of  branch  points  in  plane  zpupii ,  (j>lnd  is  given 


by 


0hid(r)  =  Im 


nf=l{x-xk)  +  i{y-yk)  1 

nfU{x-x’l)  +  i{y-Jl)\\ 


and  this  is  obviously  non-trivial  in  the  general  case.  But  since  z  <  zpupii ,  for  each  (j)n,  ^id(z)  = 
0.  It  is  immediately  obvious  that  in  the  general  case  that  each  of  the  (j>„  only  contain  least 
square  phase.  So  <j>„  <fepupil  since  <j>hid  ±  0.  Therefore,  /„  ^  /(x,y,zpupii)-  So,  an  arbitrary 
distribution  of  branch  points  cannot  be  created  in  the  propagation  from  plane  z  to  plane  z  +  £ 
with  e  arbitrarily  small. 
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5.2.  Atmospherically  Created  Branch  Points  Can  Appear  in  Pairs 

But  consider  a  very  specific  case:  suppose  that  in  plane  zpUpii  two  branch  points  with  opposite 
polarity  are  created  a  distance  d  apart.  Suppose  without  loss  of  generality  that  the  two  branch 
points  are  located  at  (0,0)  and  (0 ,d).  Let  <t>(x,y)  be  the  phase  of  a  single  branch  point.  Then 
the  phase  of  a  branch  point  pair  is 

<t>(x,y  +  r l)  —  <I>(x,y)  =  atan'— —  —  atan- 

x  x 


which  for  the  special  case  of  d  <  x  and  d  <  y 


>  0,  gives 


<f>(x,y  +  c/)  —  4>(x,y)  =  atan  ,, 

-yA  _l_  -\)Z 


dx 

x2  +y2  +  yd 


Then  it  is  trivial  to  show  that  in  this  regime 


(17) 


dx 


dx 


lim(T>(x,y  +  c/)-T>(x,y))  =  lim  ,  0  ,  „  ,  ,  , 

d- 0V  v  v  ”  d^ox2+y2+yd  3  \x2+y2  +  yd 


=  lim  d  - 


d^o  yx2+y2+yd  3  \x2+y2+yd 
0  as  d  — >  0. 


(18) 


So  pairs  of  branch  points  can  be  created  at  zpupii  such  that  @iwi  ~  0,  but  Equation  18  omits  the 
region  when  x  <  d  or  y  <  d,  i.e.  this  is  not  sufficient  to  demonstrate  uniform  propagation  but  it 
is  evocative. 

To  remedy  this,  consider  that  the  mapping  atan()  is  smooth  except  along  the  branch  cut; 
smoothness  implies  that  the  derivative  exists  and  is  finite.  So  again  consider  the  two  branch 
points.  As  is  typically  done  in  AO,  the  two  branch  points  will  be  connect  with  a  branch  cut.  The 
phase  of  the  first  branch  point  is  atan(j)  with  principal  domain  [0,  2tt] ;  this  places  the  branch 
cut  along  the  points  {(x,y}  =  ([0,°°],0).  The  phase  of  the  second  branch  point  is  atan  ( )  with 
principal  domain  \—n,n}\  this  places  the  branch  cut  along  the  points  {(x,y}  =  ([— °°,of|,0).  So 
with  atan(-)  —  atan(-7±<(),  the  branch  cuts  cancel  outside  {(x,y}  =  ([0,£/],0),  yielding  the  line 

fi(Ay;Vpii)  =  {[M]>°}- 

Since  atan()  is  smooth,  (atan(^))  exists  except  along  the  branch  cuts.  So,  as  is  given  by 
elementary  calculus, 

atan(^^"^)  —  atan(-)  =  d  atan  +  Old2)  (19) 

x  x  \xJ 


which  in-turn  implies 


atan(^  —  atan(-) 
x  x 


<dC  +  0(d2) 


(20) 


where  C  = 


lim  sup  (arctac(y))'  and  x,y  £  M./B ,  and  O(-)  is  “on  the  order  of'  Hence, 


given  e,  3D  <  —  s.t. 


W<D, 


atan('— ^ — )  —  atan(-) 
x  x 


<  ^  Vx,y  G  K2/B. 


There  may  be  question  about  choice  of  the  principal  domain.  Recall,  for  some  argument  (x,y)  £ 
R2,  atan( y)  =  0  +n2n.  In  modm  algebra,  d  +n2n  =  0,  so  choice  of  the  principal  domain  does 


not  change  0  merely  the  location  of  the  2k  discontinuity  and  the  preceeding  argument  goes 
through  unchanged.  The  restriction  to  separation  in  the  vertical  axis  is  unimportant,  and  all 
orientations  can  be  easily  obtained  via  rotation  of  coordinates. 


So,  for  a  pair  of  branch  points  of  opposite  polarity,  when  their  separation  is  small  enough. 


Zpupii) 

(21) 

\f{x,Zn)  —  /(r:  Zpupii)  | 

(22) 

/(r,z„)  —  A  (r ,  Zpupii )  el<t>lms  (r*zpupii) +'*’«<?  (r>Zpupi1 ) 

(23) 

/(r,z„)  —  A ( r , Zpupii ) el®lms (r,Zpupil )  +  A ( r , zpupii) el®lms (r,Zpupil )  _A(r,zPupii)e'<I>'mi(r’Zpupil)+4>Aw(r’j:pupil) 

(24) 


< 


/(r ,  z„)  —  A  (r , ZpUpii ) e'®,ms (r,Zpupil ^  |+|A(r,Zpupii)e'<I>,ms(r’'>Pil)  -A(r,Zpupii)e'$',“(r’Zpupil)+<I>,’'''(r’::pupil) 


< 


e 

4 


T  A(r,zpUpii) 


J  ->Z pupil) 


(25) 

(26) 


Then  examaning  the  last  term 


'^(f’^pupil) 


1  pupil) 


|A(r,zpUpii)  |  1 1  -  cos(4>/„d(r,ZpUpil))  -  i sm(Q>hid(r,  Zpupii ) )  | 

(27) 


^  4 (r. Zpupj i )  |l  cos(<I>fad(r,Zpupii))|  +  | /0  +  / sin)^,,, (r, zpupii) ) 

(28) 


then  let  d  <  D  where  D  is  given  as  above  and  D  <  <  1 
<  |A(r,ZpUpa)|  |1-1-<7C|  +  |0+(<7C)2|  (29) 

=  |A(r,zpupil)|  \dC\  +  \(dC)2\  (30) 

£  £^ 

<  — [ - <£  Md  <D  and  Vv,y  €  R2/B  (31) 

4  16  ' 


where  D  <  [ | A(r, zpUpii)  |  C]  1  £  Mr/B.  [  Recall,  that  the  zeros  of  w(x,y)  are  isolated  and  they 
exist  at  the  endpoints  of  B  ] .  Then 


fn  /(*,)>,  Zpupil)  ^ 


£ 

4 


£ 

4 


(32) 


So, 


given  £,  3D  and  N  s.t.  Vn  >  N  and  d  <D,  /„  —  /(x,y,zpupji)  <  £  \/x,y  G  R/B 


that  is 

fn  — ’  f(x,y,  Zpupn)  uniformly 

but  only  in  the  very  special  case  when  two  branch  points  of  opposite  polarity  are  separated  by 
d,  arbitrarily  small. 

That  is,  if  a  wave  propagates  through  turbulence  and  a  pair  of  branch  points  of  opposite 
polarity  spaced  infinitesimally  close  together  does  not  violate  causality. 
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5.3.  Atmospherically  Created  Branch  Points  Cannot  Appear  Singly 

Next,  consider  the  case  of  a  wave  that  is  branch  point  free  prior  to  the  pupil  but  contains  a  single 
branch  point  in  the  pupil.  Specifically,  suppose  we  have  a  wave  propagating  in  the  z  direction 
such  that,  prior  to  the  plane  zpupii,  the  wave  does  not  contain  branch  points,  and  in  plane  zpupii 
a  single  branch  point  appears.  Suppose  further,  without  loss  of  generality,  that  the  branch  point 
is  created  at  x  =  y  =  0.  Then  for  z  <  zpUpii 


(fihid  =  0  => 


and  for  z  >  zpupii 


(fihid  =  atan 


0,1 0  —  (films 

B(x,y,z  <  £pupii)  7^  ^ 


0AO  =  0/#iM  +  atan(^) 
<  Zpupil)  0 


Then,  as  before,  let 

fn  :=  /(r,z—  -)  =  A„(r)e,$"^ 
n 

with  <t>„(r)  G  (films  (r) .  As  before. 


(33) 


(34) 


* fin 


<fitms(x,y)  +  atan 


which  implies 

fn  ^  f(x,y,Zpupii). 

So  a  single  branch  point  cannot  occur  without  violating  causality. 

This  is  easily  extended  to  show  that  any  odd  number  of  branch  points  cannot  occur. 

5.4.  Arbitrary  numbers  of  Branch  Points 

Given  a  wave  with  a  creation  pair,  suppose  the  wave  propagates  further  and  additional  branch 
points  begin  to  form.  Following  the  same  type  of  proof  as  in  Sections  4  and  5,  it  can  be  shown 
that  the  second  onset  of  branch  point  formation  also  occurs  in  pairs  of  opposite  polarity  with 
infinitesimal  separation.  In  this  way,  the  pupil  can  be  fully  populated  with  branch  points. 

5.5.  Propagation  with  a  Creation  Pair 

For  completeness,  it  must  be  shown  that  once  created,  a  creation  pair  can  propagate  uniformly. 

The  proof  is  identical  to  that  in  Section  4,  with  the  exception  that  B(pc,y;z)  is  no  longer  the 
null  set,  i.e.  it  contains  the  branch  cuts  associated  with  the  creation  pairs. 

5.6.  Section  Summary 

Given  a  wave  that  is  branch-point-free  prior  to  plane  z.  A  branch  point  pair  can  be  created  at 
z  +  £  where  e  is  arbitrarily  small  as  long  as  their  separation  is  infinitely  small  and  the  branch 
points  are  of  opposite  polarity. 

We  call  this  a  creation  pair. 
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6.  Evolution  of  the  Branch  Cuts 

Suppose  a  single  creation  pair  has  formed  and  there  is  additional  propagation  with  no  further 
branch  point  formation.  Following  the  same  type  of  proof  as  in  Section  4,  it  can  be  shown  that 
fn  — >  /  uniformly  requires  that  Bn  — >  B  uniformly  where  Bn  €  {B(x,y,z,i)}  with  zn  '■=  \  as 
before.  Note,  although  B  is  changing  with  increasing  z,  this  B  is  not  a  wave  and  hence  is  not 
propagating;  we  say  here  for  discussional  purposes  that  B  is  evolving. 

So,  not  surprisingly  then,  the  branch  cuts  evolve  uniformly  as  the  wave  moves  in  z.  Uniform 
evolution  of  B  means  that  B  is  a  smooth  function  of  z  which  implies  that  branch  points  are  a 
persistent  feature  of  the  propagating  wave.  This  stands  in  direct  contradiction  to  standard  prac¬ 
tice  in  adaptive  optics  where  typically  branch  cuts  are  assigned  via  a  frame-to-frame  minimum 
distance  metric  [5],  In  short,  we  have  shown  that  the  branch  cuts  evolve  smoothly  in  z  and  hence 
in  time  since  they  travel  with  the  wave. 

Also,  see  Appendix  B  for  a  stronger  statement  on  the  condition  on  B. 

7.  Creation  Pairs,  Branch  Cuts,  and  Principal  Domains 

Each  branch  point  has  a  branch  cut.  Since  we  have  shown  that  branch  points  appear  in  pairs  of 
opposite  polarity,  we  would  expect  that  the  branch  cut  lies  along  the  branch  point  pair.  This  is 
in  keeping  with  the  observation  in  [2]  that  branch  cuts  lie  along  minima  in  intensity. 

The  principal  domain  of  each  branch  point  is  then  chosen  so  that  the  branch  cut  lies  along  the 
direction  to  the  other  branch  point.  Therefore,  since  choosing  the  line  along  which  the  branch 
cut  resides  is  equivalent  to  a  choice  of  principal  domain,  (1)  each  branch  point  in  the  creation 
pair  has  a  principal  domain  different  than  the  other  by  n,  and  (2)  each  creation  pair,  depending 
on  its  orientation,  has  a  different  base  interval  for  its  principal  domain. 

8.  Discussion 

We  have  proven  that  when  created  by  atmospheric  turbulence,  branch  points  must  be  created  in 
pairs  of  opposite  ‘polarity’  with  infinitesimal  separation.  Yet  in  real  systems,  they  are  observed 
with  physically  measurable  separation. 

Based  on  this  evidence,  three  interesting  implications  presents  themselves.  First,  the  creation 
pairs  will  have  the  velocity  of  the  turbulence  layer  that  created  them.  An  experimental  demon¬ 
stration  confirming  this  observation  will  be  presented  in  the  second  paper  in  this  series. 

Second,  since  we  have  shown  that  branch  points  must  be  created  in  pairs  infinitesimally  close 
together,  yet  they  are  observed  to  have  finite  separations  in  wavefront  sensor  data,  they  must 
drift  apart  as  the  wave  propagates.  The  forth  paper  in  this  series  will  prove  this  to  be  true  and 
also  determine  the  magnitude  of  this  drift  component. 

Thirdly,  it  appears  that  the  separation  of  the  paired  branch  points  yields  (stochastically)  the 
distance  to  the  turbulence  layer  which  caused  their  creation.  That  is,  given  a  turbulence  layer 
which  generates  an  ensemble  of  pairs  of  branch  points  some  distance  later,  the  distance  to  the 
turbulence  layer  can  be  obtained  from  the  mean  separation  between  the  paired  branch  points. 
Again  an  experimental  demonstration  will  be  presented  in  the  fourth  paper  of  this  series. 

9.  Summary 

Here  we  have  begun  a  study  on  the  aggregate  behavior  of  branch  points  and  have  shown  a  nar¬ 
row  but  crucial  result:  atmospherically  created  branch  points  are  created  infinitesimally  close 
together  in  pairs  of  opposite  polarity  and  drift  apart  as  they  propagate.  This  proves  that  atmo¬ 
spherically  created  plane  branch  points  measured  in  the  pupil  of  a  telescope  contain  information 
about  the  turbulence  layer  that  created  them.  Furthermore,  we  have  shown  that  branch  cuts  are 
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also  a  persistant  feature  of  the  wave  and  that  they  evolve  uniformly  much  as  the  wave  does 
while  propagating. 

In  short,  we  have  taken  the  first  step  in  demonstrating  that  branch  points  have  a  rich  aggregate 
behavior  previously  unmentioned  in  the  literature. 
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A.  Appendix  -  Enumertion  of  the  Derivation  demonstrating  Uniform  Propagation 

It  must  be  shown  that  the  sequence  {/„}  is  such  that  /„  — >  /(r,z)  uniformly.  To  do  so  it  must 
first  be  shown  that  f(r,z),  outside  some  radius,  R,  is  arbitrarily  close  to  zero,  i.e.  given  e  > 
0, 3R  >  0  s.t.  V|r()|  >  R,  |/(ro,z)  |  <  e.  So  consider 


/(r»z) 


(35) 


and  let 

I(r,z):=  [  /KOylNV^o-Vro. 

Ja 

It  suffices  to  show  that  V|r|  >tf,|/(r,z)|<e. 

First  note  that  ro  is  bounded  since  £2  is  compact  in  R2 .  Then,  since  r  is  unbounded. 


3R  s.t.  V|r|  >  R  s.t.  /(ro,0)  and  e'2zlr°l2  vary  slowly  w.r.t  —  | r o  -  r | 

z 


(36) 


Therefore  to  good  approximation,  both  /( ro,0)  and  e‘^r°  can  be  taken  out  of  the  integral. 
Hence, 


|/(r,z)| 


/(r0,0 WtM  [  e^'dro 
Ja 

/( r0,0)eiElr(>l2-fT  /  ei||r|pcos(0)J0 

ik  r  Jn(8(0)) 


n(p) 


<1 


/(r0,0)-  /  e(flr|pc°s(e)^0 
k  Jo.{8(p)) 


Ja(p) 


(37) 

(38) 

(39) 

(40) 


where  Q.  is  integrated  in  cylindrical  coordinates  (p.  0)  ==>  ro  •  r  =  p  |r|cos(0),  flip)  is  the 
boundary  of  £2  in  the  p  coordinate,  £l(0(p))  is  the  remaining  boundary  of  £2  over  0,  and  ]o(p) 
evaluates  the  expression  along  the  p  boundary.  Then  note  in  the  last  expression  that  since  /  has 
compact  support,  and  since  z  and  k  are  fixed,  and  since  Jnfefpj)  ^  <  +°°,  |  •  |  <  +°°  Hence, 
letting  C:=  |-|, 


I'MISo 


l/(r.,)l<4R 


given  e  >  0,3R  >  0  s.t.  V|r|  >  R,  |/(r,z)|  <  -£. 


(41) 

(42) 

(43) 
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Then,  it  suffices  to  show  that  given  e  >  0,31V  >  0  s.t.  Vr  £  M2  and  n  >  N,  |/(r,z  —  — 

/(r,z)|  <  e.  So  consider. 


(fay  _ 

\fn-f(r,z)\  =  |^-e'E|r|2/(r,z)-  H- - -e  '  7(r,z) 

lHz~},) 


!*(Z-S) 


I  /Az 


< 


Hz-}) 


1  1.1  -r'-lrpi _ ! _ 

7(r,z)-7(r,z - ) - elk"e  2  "  z(z_^/(r,z) 

n  nz 


and  using  the  triangle  inequality  and  since  0  <  z  —  1  <z - 

n 


< 


7(r,z)-7(r,z--) 


Hz-i) 

and  since  is  bounded,  EL M  s.t.  Vn  >  M, 


Ifaz) 


1  1 
n 


Hz-},) 


i 

<2£- 


Therefore, 


1, 


1 


for  z  >  1  and  Vn  >  M,  |/(r,z - )  — /(r,z)|  < 

n  A  (z  —  1 ) 


7(r,z)-7(r,z-  -) 


Now  consider  two  cases,  |r|  >  R  and  |r|  <  R  with  R  chosen  as  in  Equation  36. 
Case  1:  Irl  >  R 


|/(r,z—  -)  -/(r,z)|  <  1 

n  A  (z  —  1 ) 


7(r,z)-7(r,z--) 


1 


(44) 


(45) 


(46) 


<n§±rr+2e  (47) 

|r|  A(z-l)  2 

where  C\  and  Co  are  as  given  in  Equation  42.  Then  letting  R  :=  , 

=>  for  z  >  1  and  Vn  >  M,  given  e  >  0,EI R  s.t.  V|r|  >  R,  |/(r,z—  -)  — /(r,z)|  <  e 

Cease  I 


Case 2:  n>Mandjr|<7? 

1 


|7(r,z)-/(r,z--)| 
n 

[  /(ro,0)e'4lr°l2e‘zr°'r£/ro  —  [  /(ro.O)/2^-1/'1) 'ro  e'(z-i/")r°  rVro 
Ja  Ja 

^Vro/(ro,0)e('2™|ro|2e,™r°-r  (V2" -  1 


<  /  Vr0 


/(rO)0)|  |  ( 


Q!+j3 

£  «  —  1 


(48) 


13 


where 


a  :=i — \  7  -  |ro|2  and  /3  :  =  ikro  ■  r 


Then, 


|/(r,z)-/(r,z-- 


2  z(z  +  1/n) 


1, 


z(z+  1/n) 


(49) 


=  /  <*•() 


/(ro>0) 


1  /a  +  j3\'  1  /a  +  /3\2  1  /a  +  /3x3 


1!  V  n 


2!  V  n 


< 


3!  V  n 

3  , 


< 


^ro|/(ro,0)|^i|a  +  ^|1  +  ^|a  +  j3|2  +  ^|a  +  /3| 
/(r0,0)|  (>+/3|-l) 


n  Jn 


but  |r|  <  and  ro  €  £2  =»  a,/3  <  +°°  =>  Ja(-)  <  +°°.  So  letting  C3  :=  Jq(-)  <  +°°, 
=4-  |/(r,z)  —I(r,z—  -)|  <  -C3 

72  ft 


|/(r,z-I)-/(r,z)|<Ic3I^_  +  Ie 


(50) 

(51) 

(52) 


and  since  3N  >  2|  =>  3A  >  0,  s.t.  Vn  >  N,  |/(r,z—  4)  — /(r,z)|  <  /e  +  /e  then  letting 

P  :  =  max(N ,M) 


given  e  >  0,3f  >  0  s.t.  V|r|  <R  and  n  >  P  =>  |/(r,z - )— /(r,z)|<e  (53) 


□Cas<?2 

So,  we  have  shown  that  /„  — >  f(r,z)  uniformly  on  M.2\B  for  both  |r|  >  R  and  for  |r|  <  R.  There¬ 
fore,  free  space  Fresnel  propagation  is  uniform.  Moreover,  since  analytic  functions  are  unique 
and  given  Equation  7,  this  implies  that  both  the  amplitude  and  phase  propagate  uniformly,  i.e. 
create  sets  {A,,}  and  {<)>„}  with  A„  =A(x,y,zn)  and  <t>n  =  <t>{x,y,zn)  andz„  =  zpupii  -  4  as  before. 
Then  given  that/,,  — >  f(r,z)  uniformly  ,  A„(r)  — >  A(r)  uniformly  and  <f>„(r)  — >  O(r)  uniformly. 

B.  Appendix  -  Constraint  on  the  Evolution  of  B 

The  further  stronger  statement  can  be  made:  the  newly  created  branch  point  pair  induces  both 
real  and  imaginary  parts  to  the  wave.  Initially,  the  imaginary  part  is  zero  and  the  real  part  is 
one.  As  the  wave  propagates  past  zpupii,  the  imaginary  and  real  parts  of  the  branch  point  phase 
uniformly  evolve  away  from  zero.  This  can  be  seen  explicitly  if  d  is  taken  away  from  zero  like 

1  1 

dm  :=  —  with  in  G  Z  s.t.  —  <  D. 
m  M 
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Then  let  {<t>j '”}  (r,zpupii)}  be  a  sequence  created  using  the  dm  and  from  Equations  28  and  29, 


J  _ £^/j/d(**5£pupil) 


<  1  -COs(<I>/„y(r,ZpUpii))|  +  | i'O  +  t  sin (Q'hid (r , Zpupit) ) 


<  0  — c/C  +  iO  —  i(dC)2 

„  C 

n  fC\2 

= 

0 - 

+ 

iO  —  i  (  —  ) 

111 

\mj 

Since  C  is  a  constant,  3M  s.t.  Mm  >  M,Cjm  <  e/2 
2 

<  -  +  <  e  Mm  >  M  and  Mx,y  €  R2/B 

2  4  7 


(54) 


So, 

given  e,  3M  such  that  Mm  >  M,  1 1  —  cos(4>/,1(;(r,Zpupii))  |  <  £  Mx,y  £  K. 

and 

given  e,  3M  such  that  Mm  >  M,  |i0  — isin(<J>/,y(r,zpUpji))|  <  e  Vxj  £  R 
as  was  to  be  shown. 
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